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I.  General 

This  final  report  covers  work  carried  out  by  the  principal 
investigator  and  a  graduate  research  assistant  at  the  Department  of 
Mathematics  during  the  13  months  period  from  July  1,  1979  to 
July  31,  1980  under  Contract  F49620-79-C-0186. 

The  progress  has  resulted  in  three  technical  papers  listed  in 
Section  II,  which  will  be  submitted  for  publication  as  soon  as 
completed. 

During  the  report  period,  the  following  people  contributed  to 
the  project:  Professor  James  T.  Lo  and  Mr.  Sze-Kui  Ng. 

11*  Publications 

(1)  Transition  Probabilities  of  Homogeneous  Markov  Processes  on 
Compact  Lie  Groups,  to  appear. 

(2)  Projecting  Homogeneous  Markov  Processes  onto  Compact  Homogeneous 
Spaces,  to  appear. 

(3)  Approximation  of  Brownian  Motion  Densities  on  the  Three  Dimensional 
Rotation  Group,  to  appear. 

III.  Summary  of  Progress 

An  optimal  estimation  scheme  for  continuous -time  rotational 
processes  with  one  degree  of  freedom  was  obtained  in  [1].  The  estima¬ 
tion  problem  for  continuous- time  rotational  processes  with  three  degrees 
of  freedom  turns  out  to  be  rather  difficult.  Many  attempts  were  made 
by  the  P.I.  in  the  past  few  years  on  the  problem*  While  these  attempts 
produced  minimal  results,  many  crucial  issues  have  been  placed  in 
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perspective.  Because  of  the  theoretical  and  practical  importance 
of  the  problem,  a  concentrated  effort  was  finally  made  during  the 
report  period  to  study  these  issues. 


2. 


(1)  Transition  Probabilities  of  Homogeneous  Markov  Processes  on 


Compact  Lie  Groups: 


A  signal  process  on  SO (3) (the  3-dimensional  rotation  group)  can 


be  constructed  by  injecting  a  3-dimensional  Brownian  motion  into  ^Accession 

NTIS  G?A« 

S0(3)  [2].  It  is  known  that  such  a  signal  process  is  governed  by  DTIC  T'B 

minted 

a  bilinear  stochastic  differential  equation.  But  how  do  we  compute  JustJi'icr.t ioiu. 


the  transition  density  for  this  process?  In  the  S0(2)  case,  such  py- _ 

Distribution/ 

a  transition  density  has  an  explicit  expression  called  the  folded 


normal  density.  It  is  also  called  a  theta  function.  Then  what  is 
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a  general  folded  normal  density  or  theta  function?  This  question 


is  answered  in  a  much  more  general  context  as  follows: 


Let  {T^  }  be  a  sequence  of  finite  dimensional  irreducible  unitary 


representations  of  a  k -dimensional  compact  Lie  group  G  where  T 


<A>  -  1. 


The  Infinitesimal  operators  of  T^n^  is  defined  as 


-  fcim  ^-[t  .  -  I  ], 

1  t+0  t*'  8i^ 


i  »  1 . k 


where  {g^}  are  k  one-parameter  subgroups  of  which  the  tangent  vectors 


at  the  group  identity  e  are  orthogonal.  There  ia  a  neighborhood  of  e 


in  which  a  parameterization  x^(g),...,  x^(g)  satlaflea 


Tgn)  "  eXp(  V*eG‘ 


These  local  coordinates  (x^(g)}  can  be  extended  continuously  to  the 


3. 


entire  G  such  that  they  do  not  vanish  simultaneously  except  at  e. 

Let  denote  the  left  invariant  vector  fields  corresponding  to  g^. 
Hung  [3]  proved  that  for  a  smooth  function  f  and  the  transition 
probability  P(t,  g,  E)  of  a  homogeneous  Markov  process,  the  function 
u(t,  g)  *  I  f (a)P(t,  g,  da)  satisfies  the  Kolmogrov  backward  equation, 
with  P(0,  e,  E)  concentrates  at  e.  * 

■ . !  aijxixju<t’ *> +  !»ixi*“. *> 

1,3*=!  i— i 

(  fu(t,  gh)  -  u(t,  g)  -  l  X. (h)X  u(t,  g)ldF(h), 

■*G\{e}  i-1  1  1  J 

for  a  real  vector  [a^,...,  a^] ,  a  non-negative  definite  matrix 
k 

[a..].  .  -  ,  and  a  positive  measure  F  on  G\{e}  such  that 

ij  J- 

(  k  ? 

I  x.(g)dF(g)  <  «. 

JG\{e}i=l  1 

Our  result  is  that  the  transition  probability  P(t,  g,  E)  can  be 
calculated  by  using  the  following  fact: 

Let  Cn(t).expt|iI  =  aijA^n>Aj(n)  +  1^™ 

+ Lu)It“  - 

and  p  *  (dim.  of  T^)/(the  volume  of  G). 
n  g 

Under  certain  conditions,  there  exists  a  sequence  of  functions 


fr<t,  g,  o) 


P»Xn  4.WWe>) 


such  that  for  h6C(G), 


—  +  ■  - 
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Aim  Jh(0)fr(t,  g,  a) da  =  Jh(a)P(t,  g,  do) 

where|\r|  <1,  X*  ■  1,  Aim  Xr  =  1,  n  -  0,1,...  and  Aim  N(r)  ■  «°. 
n'— “  0  n 

r->oo  r+a> 

For  example ,  let  us  consider  G  *  SO (3)  and  the  Brownian  motion 
governed  by  the  bilinear  Ito  equation. 


dX(t) 


where  [w^. 


-dt 

-dw3 
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is  the 
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standard  Brownian  motion  in  R  .  The 

transition  density  p(t,  g,  a)  of  the  process  satisfies 


If  <*•  8.  <» 

■  jitb  M  (sin  *■ 0))  +  (^Ct'  *•  0) 

a2  a2  1 

-  2  cos  g,  a)  +  — -y(t,  g,  o)J 

p(0,  e,  a)  =  « (o) 


t  i  t 

where  g  =  (a  ,  £  ,  y  )  and  o  -  (a,  8,  Y)*  It  can  be  explicitly  written 
as 


p(t,  g,  a) 


j  2s  +  i,-iO-lH  )  .t-'f)  (co.  m.Wt'.-1- 

Sir2  ” 


where  e 


-imdpj 
rnn 


•P^  (cos  B)e 
mn 

(cos  8)e~^n^  is  the  generalized  spherical  functions  of 


order  j. 


5. 


(2)  Projecting  Homogeneous  Markov  Processes  onto  Compact  Homogeneous 
Spaces : 

2 

The  2-sphere  S  ,  the  surface  of  a  ball,  is  a  natural  state  space 

for  the  directional  processes.  It  is  isomorphic  to  the  quotient 

2 

space  S0(3)/S0(2).  A  signal  process  on  S  can  be  Induced  by  projecting 

2 

a  signal  process  on  SO (3)  onto  S  .  If  the  signal  process  on  SO (3) 

is  Brownian,  then  what  is  the  transition  density  of  the  induced 

2 

Brownian  motion  on  S  ?  This  question  is  also  answered  in  a  general 
context  as  follows: 

Let  H  be  a  closed  subgroup  of  the  Lie  group  G.  Corresponding  to 
a  projection  II :g  Hg  from  G  onto  the  homogeneous  space  HS G,  there 


exists  a  C  mapping  s  such  that  11°  s  =  1.  For  each  s,  an  element 
g£G  is  expressed  uniquely  as  s^(x)h  for  some  x€H\G  and  h6H.  A  homo¬ 
geneous  process  on  H\G  can  be  induced  by  projecting  through  II  a  homo¬ 
geneous  process  on  G. 

The  transition  probability  P(t,  x,  B)  of  the  resulting  process 
on  H\G  can  be  calculated  by  using  the  following  fact: 

Under  certain  conditions,  there  exists  a  sequence  of  functions 

£r(t'  *• 


such  that  for  h6C(H\G), 


Urn  [  h(y)f  (t ,  x,  y)dy(y) 

r-*co  Jh\G  r 


j  » 

■'H\G 


h(y)P(t,  x,  dy) 


where  y  is  an  Invariant  measure  on  H\G  and  the  other  symbols  are 
defined  as  In  (1). 
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6. 


For  example,  the  transition  density  of  the  homogeneous  process 
that  Is  the  projection  of  the  example  In  (1)  satisfies 


3p,^  x  1  3 

3?(t’  *'  y)  m7uTe  m 


x>  y) 


+  — x,  y) 
sin  0  3y 


P<0,  0,  y)  -  6(y) 


where  x  =  (a',  0')  and  y  =  (a,  0),  the  spherical  coordinates.  It 
can  be  explicitly  written  as 


P(t,  X,  y  )•  X  e’k<k  +  1)c  l  vf  -  -  Y. 


k«l 


n«-k 


k  2 


k  2 


0) 


where  “  Y»  3) »  -k  n  k  are  the  unnormalized  spherical 

functions  of  the  k-th  order. 


(3)  Approximation  of  Brownian  Motion  Densities  on  the  Three  Dimensional 
Rotation  Group: 

The  transition  densities  obtained  above  are  very  complicated 
functions,  as  illustrated  by  the  explicit  expressions  for  the  transi¬ 
tion  densities  of  the  simplest  possible  Brownian  motions  on  SO (3)  and 
2 

S  ,  In  order  to  use  them  for  estimation,  a  natural  question  is  whether 
they  can  be  approximated  closely  by  simpler  closed-form  functions. 

Let  us  restrict  our  attention  to  the  density  obtained  in  the 
example  in  (1): 


p0(a.  e»  Y) 


p(t»  ®#  (o»  3»  Y)) 


00 


1 

4-0 


24+1  -o4(4  +  1)  r  -im(a  +  y)^4  , 

Sir  e  i  •  pnn<co8 


(b)  Given  any  e  >  0  and  positive  integer  N,  there  exists  a  number 
K(e,  N)  such  that,  for 

k  >  K(e,  N),  Si2  <  Nk,  and  <  m  <_  H  , 

C*  (k)  -  +  15  <  e. 

linn  — 

It  is  noted  that  (a)  is  analogous  to  those  for  spheres  and  circles 
obtained  in  [4]  and  [5].  (b)  is  new  and  stronger  than  (a),  (b)  for 

spheres  and  circles  is  also  obtained  in  our  work. 
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